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I make a short review of the most recent determinations of the quark masses and 
run them to the Mz energy scale 


GUT and SUSY theories predict some relations among the fermion masses 
(or more properly among Yukawa couplings) at the unification scale, Mcut ~ 
10^®(GeU) . For instance, in SU{b) we have the usual lepton-bottom quarks 
unification, ht = Ht, hg = = he, or the modified Georgi-Jarlskog 

relation, hb = hr, hg = h^/Z, hd = ihe, while for 50(10) typically we get 
unification of the third family, ht = hb = hr = hi,.,.. This together with the 
RGE’s provide us a powerful tool for predicting quark masses at low energies. 

In the following I will not talk about Unification but rather try to make 
a short review of the most recent determinations of the quark masses and to 
do the running to Mz, the mass of the Z boson. Why M^?, because for 
model building purposes it is a good idea to have a reference scale, because 
extensions of the SM appear above Mz and because below Mz the strong 
coupling constant ag is really strong and then special care has to be taken on 
the running and the matching in passing a heavy quark threshold. 

Because of confinement to define what is the mass of a quark is not an easy 
task. For leptons it is clear that the physical mass is the pole of the propagator 
but for quarks we need a precise theoretical framework and the masses appear 
more like coupling constants. Dealing with quark masses we can find the 
Euclidean mass, Me{p^ = —M^), defined as the mass renormalized at the 
Euclidean point = —M^, it is gauge dependent but softly dependent on 
^QCD- However, it does not appear in the most recent works. We can talk 
about the “perturbative” pole mass, M(p^ = M^), perturbative because only 
order by order in perturbation theory the pole of the propagator is well defined. 
It is gauge and scheme independent but appears to be ambiguous because of 
non-perturbative renormalons. Finally rh{p), the mass renormalized in the 
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Table 1: Recent determinations of the light quark masses from second order QCD Sum 

Rules and lattice. 


J. Gasser 

H. Leutwyler 

XPT 

0(p4) 

'^d '^u '2,771 o Qc; N/ in—3 

iTis — m nis + m 
msfrh = 25.7 ± 2.6 

Donoghue 

et al. 

XPT 

0(p4) 

Itti oil 1 n—3 

iTis — m rris ± m 
rus/m = 31. 


Bijnens 
Prades 
de Rafael 

FESR 
Laplace SR 
(pseudo) 

NNLO 

(m„ ± md){lGeV) 

= 12.0 ±2.5 

a.(Mz) =0.117(5) 

Ioffe 
et al. 

Isospin viol, 
in QCD SR 

NNLO 

{rhd - fhu){lGeV) 

= 3.0 ± 1.0 

A = 150. 


S. Narison 

r-like SR 

NNLO 

NLO 

rhsilGeV) = 197(29) 
fhsllGeV) = 222(22) 

a,{Mz) = 0.118(6) 

M. Jamin 
M. Miinz 

QSSR 

(scalar) 

NNLO 

rhsllGeV) = 189(32) 

as{Mz) = 0.118(6) 

Chetyrkin 
et al. 

QCD SR 

NNLO 

rhsilGeV) = 171(15) 

a,(Mz) =0.117(5) 

Allton 
et al. 

quenched 

Lattice 

NLO 

rhs{2GeV) = 128(18) 

A*^ = 240. ± 90. 


Modified Minimal Substraction scheme or its corresponding Yukawa coupling 
related to it through the vev of the Higgs, fh{fi) = v{fj,)h{fi). 

For the light quarks, up, down and strange, chiral perturbation theory Era 
provide us a powerful tool for determining renormalization group invariant 
quark mass ratios. The absolute values, usually, the^mnning mass at l(G'ey) 


can be extracted from different QCD Sum Rulesl3’i0lc or for the strange quark 
mass from lattice EEl. For theJ^asw-fluarks, bottom and 

lattice calculatkmsi 


either with QCD Sum Rulesli 
soon for the bottom quark with jet physics at LEP 


e can deal 
or may be 
For the tcuL quark 

we have the recent measurements from CDF an D0 at FERMILABo I will 
identify it with the pole mass. 


I have summarized in tables and ^ all of these recent determinations of the 
quark masses. Of course the final result depends on the strong gauge coupling 
constant used in the analysis, for this reason I quote it too. In the running I will 
take a^P{Mz) = 0.118 ±0.006 for masses obtained from QCD Sum Rules but 
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Table 2: Recent determinations of the heavy quark masses from QCD Sum Rules, lattice 
and FERMILAB. 


S. Narison 

QSSR 

T 

NLO 

rht{Mt) = 4.23(4) 
fhc{Me) = 1.23+[55 
Mb = 4.62(2) 

Me = 1.42(3) 

a,(Mz) =0.118(6) 

S. Narison 

non-rel 

Laplc.SR 

NLO 

= 4.69+^ 
Me^« = 1.45+S 

as{Mz) = 0.118(6) 

Dominguez 
et al. 

rel,non-rel 
Laplc.SR 
J/'L, T 

LO 

1/ml 

Mb = 4.70(7) 

Me = 1.46(7) 

= 200 - 300 
= 100 - 200 

Titard 

Yndurain 

qq 

potential 


fhb{rhb) = 4.397;';|33| 
rhcimc) = 1.306(')|35| 

as{Mz) = 0.117(5) 

M. Neubert 

QCD SR 

NLO 

l/mq 

Mb = 4.71(7) 

Me = 1.30(12) 



Crisaf 
et al. 


i 

Lattice 
in B-meson 

mb{mb) = 4.17(6) 


Crisaf 
et al. 

A 

L6 

i 

Lattice 
in B-meson 

mb{rhb) = 4.15(7) 


Davies 

et al. 

NRQCD -1- leading 
rel and Lattice 
spacing, bb 

Mb = 5.0(2) 
rhb{Mb) = 4.0(1) 

„g^= 0.115(2) 

El-Khadra 

Mertens 

Fermilab action 
in quenched Lat 

Me = 1.5(2) 



CDF 

Mt = 176. ± S.{stat) ± 10.{sys) 

mean 

DO 

Mt = 199.+^’(;(stat) ± 22.{sys) 

Mt = 180. ± 12. 
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for lattice masses I will run with the latticed result ai^\Mz) = 0.115 ±0.002. 
This values are consistent with almost all the references. For those that differ 
an update is needed but this is beyond the goals of this paper. For instance, 
S. Narisontil makes two different determinations for the bottom and the charm 
quark masses. In the first, and for the first time, he gets directly the running 
mass avoiding then the renormalon ambiguities associated with the pole mass. 
The second one, from non-relatiidstic Laplace Sum Rules, is in fact an update 
of the work of Dominguez et al.lla. 

The 0{a‘^) strong correction to the relatitm between the perturbative pole 
mass and the running mass was calculated inEJ 


M ^ ^ 4 a,(M) ^/ a.(M) 

fh{M) 3 TT V 
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( 1 ) 


where Kt ~ 10.95 for the top quark, Ki, ~ 12.^or the bottom and Kc — 13.3 
for the charm. As pointed out by S. Narisonllil equation is consistent with 
three loops running but for two loop running we can drop the 0{a^) term. 
Recently, the electroweak correction to the relation between- the perturbative 
pole mass and the Yukawa coupling has been calculated cl. However, this 
correction is small, for instance for the top quark it is less than 0.5% in the 
SM for a mass of the Higgs lower than 600(G'eH) and at most 3.% for Mh — 
l{TeV), and for consistency one has to include it only if two loop electroweak 
running is done. 

Instead of expressing the solution of the QCD renormalization group equa¬ 
tions for the strong gauge coupling constant and the quark masses in terms of 
Age D we c^solve the running as an expansion in the strong coupling constant 
at one loopt£l. At three loops 


as(/i) = Q;i^)(/i) (l±Ci(M)ai^^(M) + C2(/i)(ai^^(M))^) , (2) 

m(/i) = + + (3) 

where (fj,) and (/i) are the one loop solutions 

(4) 

1 ± as(no)Pot 

with t = l/(47r) log^^Z/Tg, K{fi) the ratio K{fi) = and 


ci(Ai) = -h^ogKifi), 
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( 5 ) 


C2{^l) = hl\ogK{^i) [\ogK{^i) - 1] - (&? - 62 ) [1 - K{fi)], 
[ih -gi) [I-K{fj.)]+bi \ogK{g)], 

d 2 {g) = ^|[/ 3 o (&2 - bl) + 270(61 - gif] [1 - K{g)f 
+ fio{g2-higi)[l-K'^[f)\ 


+ 


47061(61 


51 ) [1 - K{f)] - 2/3o6igi +6i(/3o + 27o)logi^(^) 


\ogK{g)^, 


where 


61 


/3i 

47r/3o’ 


62 


/92 

(47r)2/3o ’ 


51 = 


7i 

7170 ' 


52 


72 

9 ’ 

7r^7o 


( 6 ) 


are the ratios of the well known beta and gamma functions in the MS scheme 


Po = 

11 

II 

Pi = 

102 - 

P 2 = 

7i = 

101 5 

12 - 18^"'’ 

72 = 


1 
2 

1249 


5033 325 , . 

2Sbl-^NF + —Nl 


277+ 180C(3),^ 35,^2 

-108- 


( 7 ) 


and C(3) = 1.2020569... is the Riemann zeta-function. Our initial condition 
for the strong coupling constant will be as{Mz)- Then we will run as from 
Mz to lower scales, i.e. for instance /zo = Mz or the upper threshold. On the 
other side, for the masses we will run from low to higher scales then we need 
the inverted version of equation |3| 


m(/ro) (l - di{g)af\g) + {dl{g) - d 2 ig))iaf\g)f'^ . 

( 8 ) 

The beta and gamma functions depend on the number of flavours Np 
therefore we have to decide ^here we have five where w^have four flavours. 
The trick is as was done bycS and recently corrected by a to built below the 
heavy quark threshold an effective theory where the heavy quark has been 
integrated out. Imposing agreement of both theories, the full and the effective, 
at low energies they wrote fj, dependent matching conditions that express the 
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parameters of the effective theory, with iV — 1 quark flavours, as a perturbative 
expansion in terms of the parameters of the full theory with N flavours 




X (fi) 1 (llx 
6 TT I2I3 2 


11 

Y 




Y) 




{^J) 


1 - 


12 


x" + 


5a; 


89 

36 


a 


N-l 


(m) 


( 9 ) 


with X = log rfi^ (/i) /, where m(/i) is the mass of the heavy quark we decouple 
at the energy scale ^ and rhi{^) are the masses of the light quarks. This 
matching conditions make the strong coupling constant and the mass of the 
light quarks discontinuous at the thresholds. However, taking the matching 
in this way we ensure, as pointed out explicitly in Ej, that the final result 
is independent of the particular matching point we choose for passing the 
threshold. As it is independent, the easiest way to implement a heavy quark 
decoupling is to take the threshold as the running mass at the running mass 
scale, i.e. jith = fhifh) or equivalently a; = 0, then the discontinuity appears 
only at two loops matching. 


Table 3: Running at the NLO and NNLO of the top quark mass to Mz, cSs'^{Mz) = 
0.118 ±0.006, ai^\Mz) = 0.117 ±0.006. 



rht{Mt) 

Tfit{rfit) 

fht{Mz) 

NLO 

172. ±12. 

173. ± 12. 

182. ±13. 

NNLO 

170. ± 12. 

171. ± 12. 

180. ± 13. 


I have summarized in tables |, I and I the result for the running of the 
quark masses to Mz- I mean by NLO connection between the perturbative 
pole mass and the running mass dropping the 0{a^) term, running to two loops 
and matching at one loop, i.e. strong gauge coupling and masses continuous 
at ^J'th = rfi{m). Three loops running and matching as expressed in equation ^ 
with a; = 0 correspond to NNLO. For consistency with the original works we 
can only do the running for the bottom and the charm quarks just to NLO. For 
the light quarks the running is consistent to NNLO using the threshold masses, 
m(m), of the bottom and charm quarks determined at NLO. I propagate the 
errors in the running in such a way we maximize them. 

It is informative to notice that the running mass of the top quark is shifted 
about 7{GeV) down from its perturbative pole mass that is of the order of its 
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Table 4: Running at the NLO of the bottom and charm quarks masses to Mz and running 
masses at the running mass scale needed for thresholds. For masses extracted from QCD 
SR af\Mz) = 0.118 ± 0.006, for lattice af\Mz) = 0.115 ± 0.002 



fh(rh) 

fh{Mz) 

S. Narison 

fhi,{fhb) = 4.29 ± 0.04 
Tnjyfhc) = 1.28 ± 0.04 

fhb{Mz) = 2.97 ±0.13 
rhc{Mz) = 0.52 ± 0.09 

S. Narison 

fh.{,(?7ib) = 4.35 ± 0.05 
Tnjyfhc) = 1.31 ± 0.06 

mb{Mz) = 3.03 ±0.13 
mc{Mz) = 0.54 ±0.10 

S. Titard 

F.J. Yndurain 

?7i6(mQ = 4.397+^:^^“ 
fhcirhc) = 1-306;°;°^^ 

fhbiMz) = 3.07 ±0.11 
ttiXMz) = 0.52 ± 0.08 

M. Neubert 

fhb{fhb) = 4.37 ± 0.09 
rhc{rhc) = 1.17 ± 0.12 

mb{Mz) = 3.04 ±0.17 
fh^Mz) = 0.45 ±0.14 

Crisafulli et al. I; 

fhb{fhb) = 4.17 ± 0.06 

fhb{Mz) = 2.93 ±0.08 

Crisafulli et al.li3 

rnbirnb) = 4.15 ± 0.07 

fhb{Mz) = 2.91 ±0.09 

El-Khadra et al. 

rhcifhc) = 1.36 ± 0.19 

mc{Mz) = 0.61 ±0.15 

Davies et al. 

rhbifhb) = 4.13 ± 0.11 

mb{Mz) = 2.89 ±0.12 

mean 

= 4.33 ± 0.06 
TTiX'^c) = 1-30 i 0.08 

fhbiMz) = 3.00 ±0.12 
fh^{Mz) = 0.52 ±0.10 


error. Therefore it is important to clarify which mass CDF and D0 are talking 
about. I have decoupled the top quark at Mz otherwise it makes no sense to 
run the top down. This fact shifts down slightly the strong coupling constant 
in Mz, from {Mz) = 0.118 ± 0.006 we get a^s\Mz) = 0.117 ± 0.006 but 
has no effect on the masses because the errors screen the difference between 
the theory with 5 and 6 flavours. Curiously the running of the top to Mz 
cancels the difference between the running and the pole mass. 

One has to be very careful in comparing the running of the masses ob¬ 
tained from QCD Sum Rules and those obtained from lattice because I took 
different values for the strong coupling constant at Mz- In addition, we have to 
remember that the error in the running is dominated by the error in the strong 
coupling constant. However imsHmpressive to notice the good agreement of 
the results obtained in lattice EIH13 with the_running of the masses from QCD 
Sum Rules. Even, the APE-CollaboratioiLl£l, has improved recently its result 
by increasing the statistics on the lattice!!^. 

We can now play a game combining the light quark masses of table || 
with the ratios obtained from xftT. The mean value of the strange quark mass 
together with the Bijnens et at.Q result gives 
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2 m, 


= 33. ±12. 


( 10 ) 


mu + md 

in agreement with the xPT result. Being conservative we can also get for the up 
and down quarks mu{lGeV) = (3. ± 2.)MeV and fhd{lGeV) = (9. ± 2.)MeV 
that translate into fhu{Mz) = {l.h±1.2)MeV and fhd{Mz) = (4.1±1.7)Mey. 


Table 5: Running of the light quark masses to Mz- For masses extracted from QCD SR 
(Mz) = 0.118 ± 0.006, for lattice (Mz) = 0.115 ± 0.002. First box is NLO running, 
second and third boxes are NNLO running. 



m{lGeV) 

m{Mz) 

S. Narison 

rhsilGeV) = 222. ±22. 

fhs{Mz) = 105. ±28. 

Allton et al. 

rhsilGeV) = 156. ± 17. 

ms{Mz) = 78. ± 15. 



m{lGeV) 

fh{Mz) 

S. Narison 

fhsilGeV) = 197. ±29. 

rhs{Mz) = 88. ±31. 

Jamin / Miinz 

msllGeV) = 189. ±32. 

fris{Mz) = 85. ±32. 

Chetyrkin et al. 

fhsllGeV) = 171. ±15. 

ms{Mz) = 75. ± 23. 

mean 

rhsilGeV) = 186. ± 30. 

rhs{Mz) = 83. ±30. 


Bijnens et al. 

{rhu ± Thd){lGeV) = 12.0. ± 2.5 

(m„ ± md){Mz) = 5.4 ± 2.2 

Ioffe et al. 

{rhd - m„)(lGeE) = 3. ± 1. 

{rhd-mu){Mz) = 1.4 ±0.7 


To summarize, the running of the quark masses to the Mz energy scale 
gives us a running top mass that is around its perturbative pole mass, mt{Mz) = 
(180.±13.)G'ey, for the bottom and the charm quark we get fhb{Mz) = (3.00± 
0.12)G'ey and fhc{Mz) = (0.52 ± Q.lQ)GeV respectively, while for the strange 
quark we have a result affected by a big error ms{Mz) = (83. ± 30.)MeV. The 
same happens for the up and down quarks, we get m„(Mz) = (1.5 ± 1.2)MeV 
and md{Mz) = (4.1 ± 1.7)MeV. 
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